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Abstract 



^SJ The Heisenberg ferromagnet has symmetry group SU(2). The prop- 

^ erty known as ferromagnetic ordering of energy levels (FOEL) states that 

^^D the minimum energy eigenvalue among eigenvectors with total spin s is 

>0 monotone decreasing as a function of s. While this property holds for 

^^ certain graphs such as open chains, in this note we demonstrate some 

. counterexamples. We consider the spin 1/2 model on rings of length 2n 

f — for n — 2,3,...,8, and show that the minimum energy among all spin sin- 

fy^ glets is less than or equal to the minimum energy among all spin triplets, 

^~^ which violates FOEL. This also shows some counterexamples to the "Al- 

^~~^ dous ordering" for the symmetric exclusion process. We also review some 

^ of the literature related to these examples. 

j^ Keywords: Heisenberg model, quantum spin systems, simple exclu- 

Vh sion process, ordering of energy levels, Aldous ordering, spectral gap, 

Temperley-Lieb algebra, Bethe ansatz. 

MCS numbers: 82B10, 81R05, 81R50. 
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1 Introduction 

We consider the spin-1/2 quantum Heisenberg ferroniagnet on even length cy- 
cles. For a positive integer N, the A^-cycle is the graph Cat with vertex set 
{1, . . . , A^} and edges 

{1,2}, {2,3}, ..., {iV-l,7V} and {l,N}. 

This is the Cayley graph for the additive cyclic group Z/(A^Z) with the generator 
set S = {1, —1}. We consider even values: TV = 2n. 
The Hilbert space for the cycle Cn is 

with single site Hilbert spaces 

Hi, ■ ■ ■ jHn — C , 

each with a prescribed orthonormal basis {|t) , \l)}. In this basis, the SU(2) 
spin operators are S^'^'^ — ^a^-^'^ for 



1 

1 



-i 
+i 



and 



+1 
-1 



Denoting the two-by-two identity matrix as 1^2 , the spin operators at the fcth 
site are 



g..,V,Z ^ ^Mk-1) ^ g^.y,, ^ ^^(„^fc) ^^^ ^ ^ ^^ 



,N . 



Then the spin-1/2 fcrromagnetic Heisenberg Hamiltonian on Cjq is 

H = hi + ■ ■ ■ + h^ , 



where for fc = 1 , . . . , A^ — 1 , 



hk = (1/4)1 - Sfc • Sfe+i 



and 



tN 



(1/4)1 -S„. Si 



We use the notation 1 for the identity operator on 'H{Cn), and the usual spin- 
matrix dot-product is 

^k ■ ^i — ^k'-'e '' ^k^i ' '-'k'-'e ■ 

We have shifted the Hamiltonian to have ground state energy equal to zero. 

The Hamiltonian operator has the symmetry group SU(2). It commutes 
with the total S*^, S*^ and S^ operators 



QX.y.z QX.y.z _. _. nx, 
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It also commutes with the total spin operator, which is sometimes called the 
Casimir operator, 

Stot = ('5'tot) + ('5'tot) + ('5'tot) • 

This is the generator of the center of the representation of SU(2); therefore, 
total spin S^^^. commutes with the "magnetization" operator 5*^^^. Since the 
self-adjoint operators S^^^ and S^^^ commute, there are simultaneous eigenspaces 

Htotis,m) = {V' e "Htot : 8^4^ = s(s + l)-!/!) and S^^^'lJJ = rmp} , 

for s G {0, 1, ... , [iV/2j} and m e {— s, — s + 1, . . . , +s}. It is here that we 
first use the fact that N is even. Writing N = 2n, we have the direct sum 
decomposition 



Htot =00 Htot(s,m) 



s— m— — s 

The maximum possible value of s is n, and similarly this is the maximum possible 
value of m. 

Defining the spin raising and lowering operators 

o± _ ex , ■ cv 
'^tot — "-"tot ^ ''-'tot I 

these also commute with H. Moreover, repeated applications of these raising 
and lowering operators may be used to map each space Utotis, rn) for a fixed 
value of s onto all the other spaces with the same s and different rn-values. 
Therefore, for a fixed value of s, the spectrum 

spec(_ff \ ntot{s,m)) 

is the same for all m. Considering an eigenvector '0 € T~l-totis^m), one often 
refers to the integer n — rn as the number of "magnons." This is the deviation 
of the magnetization eigenvalue from its maximum possible value, which is n. 
Equally important is the notion of "spin deviation," which is n — s. Let us 
denote 

Eo{C2n,k) — TmT;ispec{H \ Htotin — k,n — k)) for k — 0,...,n. 

The number Eo(C2mk) is the minimum energy on the length N = 2n cycle, 
among all eigenvectors with k spin deviates. 

We say ferromagnetic ordering of energy levels (FOEL) holds at level k € 
{0, 1, . . . , n} for the cycle C2n if 

Eo{C2n, k) < Eo{C2n,i) , for al\ee{k,..., n}. 

A famous result due to Lieb and Mattis establishes the opposite type of inequal- 
ity for antiferromagnets pT]. Lieb and Mattis's result was known as "ordering 
of energy levels," (OEL) which is why we call the ferromagnetic version FOEL. 
In [15], we proved that the FOEL property holds at all levels whenever the 
underlying graph is an open chain instead of a cycle: so the periodic edge {1, 2n} 
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is absent. This built on an earlier result of Koma and Nachtergaele, establishing 
FOEL at level A; = 1 for open chains. Also, see [T3', T^ for results about FOEL 
on open chains for higher spin models. The FOEL property for other groups 
has also been investigated by Hakobyan [S'. The FOEL property, combined 
with Bethe's ansatz, is useful for understanding the low energy spectrum of 
the Heisenberg ferromagnet, as well as the g-deformed XX Z model in one- 
dimension, as we showed in |16) . 

Moreover, FOEL at level fc = 1 is related to a conjecture due to Aldous 
[1] for the symmetric exclusion process, whose Markov generator is unitarily 
equivalent to the spin-1/2 Heisenberg Hamiltonian. Aldous's conjecture was 
recently proved by Caputo, Liggett and Richthammer [3]. More specifically, 
they proved a result which implies that FOEL holds at level k — 1 for all finite 
graphs, not just open chains or cycles. Prior to this, using a different technique, 
it was proved that FOEL holds at level A: = 1 for sufficiently large boxes [H [T2] ■ 
In a forthcoming paper, it will be proved that FOEL holds at level k, for every 
k, for sufficiently large boxes ^7\. Therefore, it is notable that we have found 
counterexamples to FOEL at level k for k = 2, . . . , 8, which is the result of the 
present paper. The method which proves FOEL at level fc > 1 only works for 
sufficiently large boxes, but this is not just a deficiency of the proof technique 
because in fact FOEL at level A: > 1 does not hold for all finite graphs. This also 
gives some counterexamples to the more general property of "Aldous ordering," 
introduced by Alon and Kozma in 2]. 

Main Results: For each value of n ~ 2,3, ... ,8, 

Eo{C2n,n) < Eo{C2n,n~l), 

and there is .strict inequality for n > 2. 

For n — 2 and 3, we prove this result by algebraic diagonalization of the 
Hamiltonian operator. That occupies Section 2. In Section 3, we extend these 
results using exact numerical diagonalization, which establishes the result for 
n = 4, . . . , 7. In Section 4, we discuss the relation to the Bethe ansatz, and also 
how to infer the result for n = 8 from previous results of Dhar and Shastry [7] . 

2 Exact diagonalization of the graphical model 

For doing calculations, we find it convenient to use a graphical basis, related 
to the Temperley-Lieb algebra. At first, we will consider the model on an open 
chain with 2n sites. Later we consider how the periodic boundary conditions, 
and especially the edge {1, 2n} alter the set-up. 

Starting from the all |t) spin state, one may make a fc-spin deviate vector 
as follows. For each ordered set of 2k distinct sites (ai, bi, . . . ,ak,bk) define 



*(ai,5i,...,afe,6fe) = H '^a.-.bjtr) ' 
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where Iff) is the all-up-spin state, and 



a,b 



'^a ^b ■ 



We can write this instead as 



*(ai,6i,...,afc,fefc) 



iV'av.f),- fXl 



It). 



i=i 



l^{ai,bi,...,ak,bk} 



where ipa,b — \i)a ^ \t)b — It) a <8) \i)b- There are some graphical rules, which 
are designed to obtain the right number of vectors this way to form a basis for 
'Htotin- k,n-- k): 

• For each j, aj < bj. 

• Each j G {!,..., 2n} must appear in the list (ai, 61, . . . , a^, 6j) at most 
once. If it does not occur at all it is called an "unpaired site." 

• Drawing arcs from aj to bj for j = 1, . . . ,k, no two arcs may cross. 

• No arc may span an unpaired site. 

Graphically, we represent this as follows. We draw a line of 2n dots. Above the 
line, we draw an oriented arc from aj to bj for each j, and we either draw the 
t vectors above each unpaired spin site or else we draw a (J) above the line, 
indicating the ^ spins, and draw strands from each unpaired site to this: 



*(2,3;5,8;6,7) 






More generally, for an unpaired site, we could place an f or a ,|, spin. However, 
we get a basis vector in "Htotin — k,n— k) only if all unpaired spins are 'f spins. 
Such a vector is called a "highest weight" vector because for total spin s ~ n — k 
it has the highest possible value of the S^^^ eigenvalue m, namely n -- k. 

Co- vectors, i.e., linear functionals, have a slightly different form. We define 

i^\a,b) = a(tl®6a|-6(tl®aa|, 

which has the opposite orientation so that. 

This is convenient in the graphical representation for the following reason. We 
denote tpl^i^ as an arc between the sites a and b. We also denote the covector (t| 
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or (4,| at a site by placing the spin configuration under a site. Then we have the 
useful formulas: 



= 5„ 



a 



and 



These formulas may be checked by expanding each singlet vector ipab and covec- 
tor ■0^^. The graphical interpretation is that one may contract the middle arcs 
to simplify the diagram. This is the first instance of a general procedure related 
to the graphical representation, interpreting algebraic relations as topological 
simplifications. This is also the idea behind the Temperley-Lieb algebra. 



2.1 Temperley-Lieb algebra 



We know that ^g^ij4'ab = ^2. This is represented diagrammatically as saying 
that the "loop" evaluates to —2: 



-2. 



This is the scalar inner-product. But we may also take the outer product to ob- 
tain an operator or endomorphism. The Temperley-Lieb generators are defined 

as Ui,. . . ,U2n-i, where 



Uk = lc2 «) (V-fc.fc + l^^ + l) 



, 2n-l-fe 



One may check that 



Uk = -2hk , 



for each fcG{l,...,2n— 1}. In order to represent these graphically, we represent 
a factor of lc2 as a straight vertical line segment, with the interpretation that 
contracting the line segment does implement the Kronccker delta, as it should: 



Uk 



k fc + 1 



fc + 2 



2n 
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for k — 1, . . . ,2n— 1. The Temper ley-Lieb algebra is the algebra obtained from 
these generators, using the relations 

Ui = -2Uk 
UkUk+iUk ~ Uk 



UkUk-iUk = Uk 
\k-£\>l => UkUt = U,Uk 

With the representation given above, where Uk = —2hk, these relations are all 
satisfied. The graphical interpretations of these relations are as follows, where 
we leave out unnecessary tensor factors of the identity to remove extraneous 
vertical lines: 



o 



Ky 



Ky 



r\ 






and 









r^ 



hy -uh 



These may also be interpreted as topological simplifications. The model for 
many of these algebraic formalisms of topological simplifications is given by 
Reidemeister's "moves" for knot and braid diagrams. Indeed, the Temperley- 
Lieb algebra played an important role in the discovery by Vaughn Jones of 
his knot polynomial. We recommend the article of Frenkel and Khovanov for 
the description in the setting of the quantum group SUg(2) [5]. The graphical 
representation has been important in proving FOEL for open chains, as in |15) 
and more recently |13l. 

Since Uk — — 2ft,fc, for k — 1, . . . , 2n— 1, this gives a convenient representation 
for the Hamiltonian. One may consider the formal vector space spanned by 
the diagrams for vectors, and use an operator A representing the sum of the 
Temperley-Lieb generators. Then, after finding the eigenvalues and eigenvectors 
in this basis, one may transform back into the standard (Ising) tensor basis of 
It) and 14,) vectors by expanding each t/jab singlet vector. As some examples of 
the diagrammatic rules for the graphical representations of Uk , we consider the 
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following 



U2 I ^12 ® (g) It), I = 

j=3A 



U2 U: 



23 



It), 



J = l,4 




= ^23 ® (g) It), , 
, = 1,4 



-2^-23 ® 



It) 



, ' 



, = 1,4 



U2 (-012 «) -034) 




V'M «) 0: 



23 : 



Also, since the singlet is antisymmetric, while |t) It) is clearly symmetric, this 
implies 



C/l(lt)l®lt)2) = 




= 0. 



A more general identity of this sort involves the Jones- Wenzl projector, which 
is a symmetrizcr for two or more spins. We will not need the more general 
identity. 

2.2 Graphical representation on the ring 

We want to consider the graphical representation on the ring. However, momen- 
tarily let us continue to consider the spin-1/2 Heisenberg model on the chain, 
so that the edge {1, 2n} is still absent, as is the interaction ft,2n- 

For the action on the k spin deviate subspace, we may define another vector 
space V which is the set of formal linear combinations of the arc-diagrams on the 
chain. Then we may define an operator A acting on V ^ obeying the appropriate 
rules one infers from the Temperley-Lieb algebra. This will be the sum of all 
Temperly-Lieb generators: 



A = C/1,2 + 



U' 



2n-l,2n 



(We now write Uk as Uk.k+i to emphasize both vertices involved in the arcs. 
Finally we may define a linear transformation between vector spaces. 



L -.V -^ Htot {n~ k,n— k) . 
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which implements the original definition of the arc diagrams in terms of linear 
combinations of Ising spin configuration basis vectors. Because the Temperley- 
Lieb algebra is related to the Hcisenberg model by the identities Uk,k+i = —'^hk, 
for k — I, . . . ,2n — 1, we then have the intertwining relation 

LA = -2HL . 

Therefore, if there is some eigenvector (j) of A, such that A(p ~ X(f>, then the 
vector tp = L(j> must satisfy Hip — —{l/2)\ip. This shows that —(1/2) A is an 
eigenvector of H, as long as ip = L(j) is not zero. For the case we are considering, 
L is actually an isomorphism, as can be seen by counting the number of basis 
vectors. See, for example, [TS] or [T3] for a proof of the combinatorial identity 
showing that the dimensions of V and Htot{n — k,n ~ k) match. 

So it never happens that L(j) = if (/> is nonzero. Therefore, the entire 
spectrum of {H \ T-Ltot{n — k,n — k)) can be obtained merely by multiplying all 
points of spec(A) by —1/2. 

In order to consider the cycle we now need to include /i2n- This is not a simple 
Temperley-Lieb algebra generator. We note for example, that the singlet vector 
satisfies 

V'2n,l = -(^12 + ?A23 H h V'2n--l,2n) , 

with a similar formula for the covector. But in the graphical model, we may 
merely add another generator which acts in the correct way, according to the 
rules of the Templerley-Lieb algebra. 

In other words, instead of decomposing an arc such as that representing 
'02n,i as a linear combination of the other more basic arc diagrams, we merely 
append that arc as part of the vector space V. Thus we enlarge the vector space 
y, and therefore the intertwining operator L : V^ — > T-Ltotin — k,n — k) will no 
longer be an isomorphism. It will map onto 'Htoti''^ — k,n ^ k), but L will have 
a kernel. Taking this into account, and defining A = U12 + . . . U2n-i.2n + C^2n,i, 
we have 

spec{H \ ntot{n-k, n-k)) = {-(1/2)A : 3(j) e V , Acj) = Xcj) and cp ^ ker(i)} . 

When working on the cycle, we must specify the orientation of each arc, 
because it is not uniquely determined as left-to-right as for the chain. Reversing 
the orientation of a single arc corresponds to multiplying a vector m'H{n—k,n — 
fc) by — 1, after the action of L. Note, however, that reversing two arcs does not 
change the sign. For this reason, in a Temperley-Lieb generator Uk.k+i reversing 
the orientation of both arcs, simultaneously, does not change its action. 

In principle we could allow V to be the formal linear combination of all 
oriented arc diagrams. However, we will typically choose a subset. For each arc 
configuration, we will want to consider all translates of it, under the rotation 
group which is a symmetry of the operator A (and H). This simplifies the 
calculations. However, we often times do not do the same for the reflections, 
since the full dihedral symmetry I?2n is not much more useful than the cyclic 
symmetry C2„. 
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2.3 The square C4 

We begin by considering C4 and k — 1. 

2.3.1 C4, 1 spin deviate 

We consider V to be the formal span of the following diagrams. 




!^23 




^34 




012 




" 1 " 


023 




T 


034 




r2 


_041_ 




y3 



and 



We enumerate the vertices from 1 to 4 starting at the top right, and proceeding 
clockwise. We let T be the rotation one unit to the right, so that 



»12 ■ 



Then, letting U12, C/23, C/34 and Un be the four generators in the extended 
Temperley-Lieb algebra on the ring, we have 



^12 ■ 



Therefore, for 

A = U12 + C/23 + C/34 + U41 , 

we have A(j>i2 = {—2 + T + T^)(pi2. We define translation eigenvectors next. For 
z € {l,i, —1, ^i}, we define 



^12' 




'-2 


U23 
t/34 


012 = 


T 



c/41 




T\ 



z)^Y. ^'^* 



fc=0 



These are also eigenvectors of A: 

A0(z) = (-2+z+z3)0(z) =^ 



yl 



-2 



-4 



\ r 0(1) ■ 
m 

/ L0(-i). 



(2.1) 
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Next we must calculate the transformation to the 1-spin deviate subspace 
for C4, 

We denote one-magnon Ising basis vectors in Htot for C4 as 



l-i 



and 




$„ = 




v[/„ = 




1-4 



We let T also denote the implementation of the same rotation operation on 
?^tot, through its Ising basis. Then 



*i 



Also, 

L(j,i2 = *i-*2 = (1-T)*i. (2.2) 

We define translation eigenvectors in "Htot- For z e {1, i, —1, — i}, we define 

3 



r*ii 




" 1 " 


*^ 




T 


*8 




y2 


*4. 




y3 



^(z) = ^z'^r'^^i. 



fe=0 



These are linearly independent vectors. Moreover, from (2.2 1, we have 

L${z) = {l~z^)^Jj{z). 
Therefore the kernel of L is one dimensional, spanned just by <^(1). So, using 



(|2.1|), we see that the spectrum of H on the 1 spin deviate subspace for C4 is 
spec(2i7rHtot(l,l)) = {2,4}. 



given by 

In particular. 



Eo{Ci,l) = 1, 



(2.3) 
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2.3.2 C4, 2 spin deviate 

We consider the formal vector space V, spanned by two vectors 



^12, 34 — 



Then 



Also, 




-'14,23 



© 





<p 


12,34 




1 


/A 




_014,23 




yjr 


U12' 






-2 


t/23 


012,34 — 




T 


c/34 




-2 


U41 








T 



^12,34 • 



^12,34 ■ 



We define the translation eigenvectors, for z E {1, — 1}, 



0(Z) = ^z'^T'^012,34. 
fe=0 

Then for A = U12 + U23 + C/34 + U41, we have 



A(l){z) = (-4+2z)0(z) 



A 



"-2 

-6 


) 


■ 0(1) ■ 
.0(-l). 


= 


'0 




(2.4) 



We next need to calculate L. The 2-niagnon subspace of Htot, consists of 6 
basis vectors: 



*,fc = S7S 



-]k 



j '-'k 



for 



1 < j < fc < 4 , 



where Iff) is the all-up-spin vector. We define two representative Ising basis 
vectors from "Htot, one for each translation orbit: 



^ 



12 





and v]/j3 = 
Then, using the same symbol T for translations in "Htot, 

i012,34 = *13 - *23 - *14 + *24 

= *13 - T*i2 - T3*i2 + T*i3 

= -(T + r3)*i2 + (l+T)*13- 



(2.5) 



We define translation eigenvectors. But there arc now two different sizes for the 
orbits. 
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• For z E {1, i, —1, —i}, we define 

3 



fe=0 

For z E {1, —1}, we define 



^2(2) = ^z^'T*'-* 



13 ■ 

A;=0 



These are all linearly independent vectors in "Htot for the cycle C4. From (2.51 
we see that 

L0(Z) = -(z + z3)V>l(z) + (l+z3)^2(;2), 



for z E {—1, !}■ This has trivial kernel. Therefore, from (2.4) we may read off 



the spectrum for H restricted to the 2 spin deviate subspace for C4 

spec(2i? rHtot(0,0)) = {2,6}. 

In particular, 

EoiCi,2) = 1. 



Comparing to (2.3) we see that we have the "accidental degeneracy," 

EoiCi,2) = EoiCi,l). 

This satisfies both FOEL at level 1, because £'o(C'4,l) < £'o(C'4,2), and our 
"opposite surmise," i?o(C2n,n) < EQ{C2n,'n— 1), because of the equality. This 
is the first example from our main results. 

2.3.3 Plots for C4 

In Figure 111 we have numerically diagonalized 2H, using the Ising basis in 
Matlab. Matlab uses a standard implementation of the exact numerical diago- 
nalization algorithm "Lanczos," which is considered to be reliable for matrices 
of the sizes we considered. We consider the square which is Cn for TV = 4 in 
that figure. All of the eigenvalues we have calculated here are displayed there, 
and the "accidental degeneracy" is circled. 

2.4 The hexagon Cq 

The Hubbard model on the hexagon was investigated by Heilmann and Lieb with 
a view towards degeneracies [TO] . Since the Hubbard model is more complicated 
than the spin- 1/2 Heisenberg model, this suggests that the investigation of the 
Heisenberg model may also yield interesting results. As a remark, we mention 
that they found many "accidental" degeneracies, i.e, degeneracies which are 
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not obviously related to a symmetry of the Hamiltonian. It was this which 
motivated the later work ^^. 

Let us mention the easiest results. We know that £'0(6*6, 0) = 0, because 
that is the ground state energy. That is also the only eigenvalue of H in the 
one-dimensional subspace "HtotlS, 3) of 0-spin deviates. 

The spectrum of H in the 1-spin deviate subspace 'HtotlS, 2) is also easy. We 



may work in analogy with Subsection 2.3.1 We obtain for Cq, that there are A 
eigenvectors, ${z) e F for z e {l,e'''/^ e2'''/^ -1, e^''*/^^ gSWS}^ 



z) =Y^ z^T'' 




fe=0 



and 

A^{z) = {-2 + z + z^)^{z) . 

However the kernel of the intertwining operator L : V ^ Htot(2, 2) consists of 
0(1). Therefore, for Cg, 

spec(2iJrHtot(2,2)) = {1,3,4}. 

More specifically, 

(2iJ-l)i0(e±'^*/^) = (2iJ - 3)i0(e±2^'/^) = {-2H - 4)L^{-1) = 0. 

In particular, £'o(C6;l) = 1/2. Another interesting fact relates to the energy 
eigenvalue 2H(f>{—l) — 40(— 1). As we will see, that eigenvalue of 2H also has 
an "accidental degeneracy." 



2.4.1 Cg, 2 spin deviate 

We now consider the case s = 1. This is the 2 spin deviate subspace for Cg. 
The diagrams we consider have 2 arcs. We define three basic diagrams 




^12,34 




and 



'12,45 
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With T denoting the fundamental clockwise rotation, we may write 

-2 



A 



^14, 23 
^12,34 
^12,45 



2 + T + T^ ■ 

1 -4 1 + T^ 

1 + T + T^ + T-* -4 



^14,23 
^2,34 
^2,45 



(2.6) 



We define translation eigenvectors based on these three translation orbit repre- 
sentatives. The orbits have two different sizes. 

• For z e {l,e"/3^e2W3^_i^e4W3^g5W3}^ .^^e define 

5 5 

M^) = Y^z'^T'^cj^uaz and 02(^) = ^ z'=T'=</.i2, 

• For z e {l,e2^*/3,e''''*/3}, we define 

2 
03 W = E^'^'^12. 



34 



/c=0 



;,45 



fc=0 



Then (2.61 may be used to deduce the action of A in each translation eigenspace. 



For z e {e'"*/^, -1, e^"/^}, we have 
[xi 2:2] A 



02(4 



Xi X2\ 



-2 2 + z + z'"^' 
1 -4 



'0i(^)' 

02(4. 



• For z e {1, e'^'^*/'^,e"''^*/"^}, the subgroup determined by the equation z^ 



[xi X2 X3] A 



01 (^y 

02 (^) 
.03(4 



= [xi X2 X3] 



-2 2 + z + z2 

1 -4 2 + 2z 

1 + z2 -4 



■01 (^)' 

02(4 
.03(4. 



To diagonalize A in each translation eigenspace (sector) one must solve for the 
left eigenvectors [xi,X2] or [xi,X2,X3\ of the matrices displayed above, on the 
right-hand-sides of the equations. From this we may deduce the spectrum of A 

• For z G {e'^'/'^, — 1, e^'^'/^}, which means z^ — —1, we have: 

^ - (-3 ± Vs + z-z^^^ [(^i(z) + (-1 ± V3 + 2-22) 02(z) 



= 0. 



For z e {e'^"/^, e'^^i/Sj ^^ i^ave: 

[A-i-3)]\M^)-4>2{z) + 2z^3{z) 
-7±vT7^ 



A- 



0, and 
3 ± Vl7) (/.i(z) + 4^2(2) + (1 T VT7)Mz) 



('3±\/l7] 



= 0. 
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For z = 1, we have 



Mz) + 2Mz)+2(f>3iz) 



and 



A 



(-5 ± Vtj^ [01 (z) + (^-3±Vtj Mz) + 2 (l T ^) Mz) 



Next we compute L in this basis. We define some basic Ising configuration 
vectors 



*12 



t i 



t ; 



and 



^1 






Using T for the usual shift, one may calculate the action of the intertwining 
operator L : V ^ Htot(l, 1) for Cg, 



^14,23 
^12,34 
^12,45 



l + T^ -(l + T) 

-T l + T -1 

-{T + T^) l + T 



*12 

*13 
*14 



(2.7) 



We define translation eigenvectors, as follows. 

• For z e {1, e^*/3^ g27r»/3^ _;l^ e^''*/^^ e^""'/^}, we define 

5 5 

7^i(z) = ^z'^T'^vI/12 and ^2^ = ^ z'=r'=^I'i3 . 



k=0 



fc=0 



For z e {l,e2"/3^e''"/3}, we define 



^3(^) = ^z'^T'^vI/,4, 



fc=0 



Then, from (2.7 1, we have the following formulas for L in each translation 
eigenspace. 



For z G {e'^*/'^,— 1,6^^^'/^}, which means z'^ = — 1, 



[xi X2] L 



Mz). 



[xi X2] 



l-z -l + z^' 

Z^ 1 - z2 



Mz). 
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For z G {1, e^TTi/a g47ri/3|^ which means z^ = 1, we have 



[Xi X2 X3] L 



k{z). 



= [Xi X2 2:3] 



V + z 


-l-z2 







"^l(^) 


2 

Z 


l + z2 


-2 




V'2(z) 





-z2 


l + z2 




Lv'3(^) 



From this, we can see that the kernel of L is spanned by the vectors: 

(^i(l) + 202(1) + 203(1) , <^i(-l) - 202(-l) , 0i(e±"*/3) + 02(6*^^^/3) ^^d 
0^(e±2-/3) _ 02(e±'"^/') + 2e±4'''/3^3(e±2-/3) . 

Each of these vectors is an eigenvector of A. The eigenvalues associated to 
these eigenvectors, which need to be removed are 0, —1, —3, —4. This leaves 
the following eigenvalues of H in the 2 spin deviate sector for Cg: 



a(2if fHtot(l,l)) 
In particular, 




, 2, 5-^5, 5, 5 + ^5, - 



17 



^o(C6,2) = 



7- V17 



2.4.2 Cg, 3 spin deviate 

We now consider singlets. Graphically, this corresponds to total matchings. We 
consider two representatives of translation orbits: 



'14,23,56 




and 



012,34,56 



o 



Note that for 0i4,23,56j rotation by it changes the orientation of the arc con- 
necting 1 and 4. In principle, we should also introduce the basis vector 02i,43,56 
which is related to 0i2,34,56 by reversing all orientations. However, we find that 
the translation on 0i2,34,56 does not reverse orientation. Moreover, we are not 
trying to exploit a "good signs" condition, using the Perron-Frobenius theorem. 
Therefore, we merely allow for the possibility of taking — 0i2,34,56 instead of 

?^12, 34,56- 

Letting T be the usual translation. 



''14,23,56 
^12,34,56 



1 



-4 



rp4: 



2-2T 
-6 



914,23,56 
^12,34,56 



Note the generalized matrix entry 2 — 2T, rather than 2 + 2T. This "bad sign" 
in an off-diagonal entry is precisely because we did not also include the vector 
021,43,56 which reverses all three orientations from 0i2.34.567 but rather simply 
took this to be — 0i2, 34,56 • 

We define translation eigenvectors. 
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For z G {1, e"/^, e^"/^, -1, e"^""'/^, e^""'/^}, wc define 



hiz) = ^z'^T^ 



14.23.56 ■ 



fe=0 



For z e {1, —1}, we define 



kiz) - ;^ z'^TVi2,34,56 . 



fe=0 

Then we have the following action of A. 

• For Z e {eW3^g2W3^g4^j/3^g5W3}^ 

^01 (z) = -#i(z), 

• For zG {1,-1}, 



[xi X2] A 



'M^) 

M^). 



= [xi X2] 



-4 2(l-z)(l + z2 + z4y 
1 -6 






From this we deduce the spectrum of A, as follows. 

• For Z e {eW3^g2W3^g47rj/3^g57ri/3}^ 

[A-i-4)]M^) = 0. 

• For z = — 1, 



A- -5± V13 



[(/.l(z)+(^-l±Vl3J(/.2(z) 

• For z = 1, 

[A - (-4)] Mz) ^ and [A - (-6)] [0i(z) - 2$2{z) 

Next we calculate the action of L. We define the Ising basis vectors 
t i t -I- 

1-123 = t/ U , *124 = t/ n 



0. 



= 0. 



t i t i 

*125 = l4 n and *i35 = ;i it 
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Then 



L 



'14,23,56 
^12,34,56 



T^2 T^5 T^2 I rjib 1 T^3 1-1-7^ 

-T - T-^ - T-"^ T + T^ + T"' l-T 



'*123' 
*124 
*125 
*135 



We define translation eigenvectors as follows. 

• For z e {1, e"^/'^, g^'^^/s^ _i^ e'^'"^/'^, e^'"'/'^}, we define 



V^2(z) 



= ^ z'^r'^ 



fe=0 



*123 
*124 
*125 



For z e {1, —1}, we define 

1 



135 ■ 



fe=0 



Then we may rewrite the action of L as follows. 
• For z e {eW3^g2W3^g4W3^e5''*/3}^ 



Mz) 
Mz). 



For ze {1,-1}, 



[xi X2] L 



'Mz) 
Mz) 



= [(1 - z)xi X2] 



—z z 1 —3 

— z z 1 — z 



'V'lC^y 

V^2(Z) 

.^'4(2). 



From this, we may deduce that the kernel of L is spanned by the vectors 

^(1) and (^(e^^'^^/S) . 

Removing these vectors as eigenvectors of A leaves the spectrum of H in the 3 
spin deviate space for C%: 

spec(2i7- \ Htot(0, 0)) = {5-\/l3,4,6,5 + \/l3} . 

In particular this means that 

-£^o(C6,3) = . 
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Comparing, we see that 

£;o(C6,3) = ~ = 0.697224362 < £;o(C6,2) = ~^ = 0.719223593. 

So this gives a violation of FOEL at level 2 since n = 6 implies that spin s = 1 
is {n/2) — 2. This is the second example of our main results. 

2.4.3 Plots for Ce 

In Figure IT] we have plotted the results of having numerically diagonalized 2_ff, 
on the ring Cn for several values of TV. The hexagon is TV = 6. All the eigenval- 
ues we mentioned are represented there. In addition to the violation of FOEL, 
there is also an "accidental degeneracy." The eigenvalue 4 occurs both in the 1 
spin deviate subspace 'Htot(2, 2) and the 3 spin deviate subspace Htot(0, 0). The 
violation of FOEL is indicated with an arrow, and the "accidental degeneracy" 
is circled. 



3 Counterexamples through exact Lanczos di- 
agonahzation 

In this section we report on results of numerical diagonalization of the spin-1/2 
Heisenberg model on cycles Cjv for various values of N, both even and odd. 
We used Matlab's built-in "eigs" command which is an implementation of the 
Lanczos iteration scheme. For iV = 4, . . . , 8, we calculated the full spectrum. 
This is plotted in Figure [l] For N = 10, 12, 14, we have just plotted the lowest 
eigenvalues, and zoomed-in to see the relative ordering of eigenvalues. This is 
plotted in Figure [2j For TV = 15, we have plotted the eigenvalues separately, 
since size limitations required us to plot only a subset of the relevant information. 
This is shown in Figure [31 We now explain, point by point, the various features 
of relevance. 

3.1 Details of the implementation 

For Cat, we define n = [-/V/2J . When constructing the matrix for H in Matlab, 
we restricted to the n-magnon subspace. This is the subspace 

nZ!'{{N/2) -n) = {^^ Htot : ^t^ot^ = [{l/N) - n]^} . 

Since total spin also commutes with H and S^^^, it is useful to consider the 
further direct sum decomposition: 

n 

'HZ'{{N/2) -n) = 0Htot((A^/2) - k, {N/2) - n) . 

k=0 

Since spec(if \ Htot(s, m) is actually independent of m, we see that in T-L"ot^{{N/2)- 
n) we have representative eigenvectors corresponding to every eigenvalue of H 
on the full Hilbert space "Htot- So this is the optimal 5'^ eigenspace to consider. 
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12 3 
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7V = 8 
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7- O 
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5 - G 



- 
12 3 

N = 7 



Figure 1: (Color online.) We have plotted spec(2_ff \ Htot{s,s)) for the cycles 
Cn- for N = 4, 5,..., 8. The horizontal axis is s for s G {(L/2) — k : k = 
0, . . . ,n} where we define n = [N/2\ . Each plot point is a circle with a sector 
inscribed, which corresponds to the translation eigenvalue(s) of the associated 
eigenvector(s). The dashed circles represent "accidental degeneracies," while 
the arrows indicate violations of FOEL. 
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G 
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G 
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o 
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G 
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3 4 5 6 7 

A^ = 14 



Figure 2: (Color online.) We have plotted spec(27f \ 'Htot{s,s)) for the cycles 
Cn for N — 10, 12, 14. The horizontal axis is s for s G {n — fc : k — 0, . . .n} 
where we define n — N/2. Each plot point is a circle with a sector inscribed, 
which corresponds to the translation eigenvalue(s) of the associated eigenvec- 
tor(s). In order to compare eigenvalues, we have put a dotted line emanating 
from each eigenvalue. By zooming in, one can see that FOEL is violated, as 
claimed in the main results, for these examples. 
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Energy Vs TranElalion For L=15 




-0.4 -0.2 0.2 

TranElaHon Eiganvalues 



Figure 3: We have plotted spec(2iJ \ HtotisTs)) for the cycle C15. However, 
now the horizontal axis denotes cos 61 for the T eigenvalue e*^. In this plot, we 
do not plot the total spin. But according to Sutherland's surmise, the minimum 
energy among all eigenvectors with Tip = e±27rifc/i5 jg eqya,l to the minimal 
energy among all eigenvectors with total spin equal to (15/2) — k, for k = 
0, 1, . . . , 7. Following this conjecture, we see that FOEL is satisfied for N = 15. 
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We used Matlab to do exact diagonalization. We also calculated the total 
spin of the eigenvectors. Note that the translation operator T also commutes 
with iJ, S'^Qj and S^^^. So we also calculate the translation eigenvalues e*^ for 
each eigenvector of H. All of this is plotted in Figures 1 and 2. On the horizontal 
axis we plotted total spin s G {{N/2) — k : k G {0,1, . . . ,n}, such that for an 
eigenvector ip oi H, we have 

On the vertical axis we have plotted the eigenvalues E such that 

2Hi: = Ei} . 

We plotted this in order to better compare to the A eigenvalues from the last 
subsection. It also matches the convention used by physicists, which we will 
review in the next section. 

The plot points are open circles, and in each circle we have plotted rays 
associated to all the unit complex numbers e*^ such that T^p = e*^^'- For 1 and 
— 1 the eigenvectors are typically non-degenerate. But since H is translation 
invariant and real, if < < tt and there is an eigenvector ■0+ such that 

2HtP+ = E^+ and TiP+ = e+*V+ , 

then taking ■(/;_ to be the complex conjugate of tl)-^-, we have 

2H^^ = Eij)^ and T^^ = e"'V_ . 

For N — 15, we were more limited in our use of Matlab. Therefore, we did 
not plot the total spin. We merely plotted the energy versus the translation 
eigenvalue. In the next subsection we will explain why this also gives useful 
information. 

3.2 Some conjectures 

According to our numerical examples, we seem to see the following pattern: 

• For even length chains Cjq with N = 2n, 

Eo{C2n,n) < Eo{C2n,n^l). 

We found this to be true for n = 2, ... ,7 , and as we will describe, one 
may also infer this violation of FOEL also for n ~ S. 

• For odd length spin rings, N = 2n-|-l, we did not find any violations of the 
FOEL ordering. Specifically, we found EQ{C2n+i, 1/2) < -E^o(C2n+i,3/2), 
which verifies FOEL. This might simply be that we have not considered 
sufficiently long chains. 
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In order to explain the A'' = 15 example, let us mention the following. Physi- 
cists have considered a problem related to these counterexamples, although not 
in the context of proving or disproving FOEL. In particular, for the ring Cat, 
Sutherland [TB] noted the following trend. Note that the T eigenvalues are of 
the form e^^ for 9 = ±2Trk/N for A: = 0, 1, . . . , n, where n = [iV/2j . Let 

for k G {0,1, . . . , n}. Then Sutherland's observation is 

minspec(i7 \ nl'J'^"' {k)) = minspec(ff \ Htot{{N/2) ~ k, {N/2) - k)) , 

for k e {0, 1, . . . ,n}. To the best of our knowledge, Sutherland did not prove 
this conjecture, and we also have not obtained a rigorous proof. But Suther- 
land's surmise is verified in all our examples A^ = 4, 5, 6, 7, 8, 10, 12, 14 as direct 
inspection of Figure [l] and Figure [2] reveal. 

In Figure [3] on the horizontal axis we plotted cos 9 for the T eigenvalues e'^. 
Therefore, the lowest band of energy eigenvalues in Figure |3] should correspond 
to -Eo(Ci5, k) for k = 0, 1, . . . , 7 reading from right-to-left (for increasing k in 
cos(27rA;/15)). Therefore, this does not appear to violate FOEL. 

Finally, Dhar and Shastry reconsidered some of Sutherland's results in [7]. 
In their Figure 2, they considered the spin ring Cig, and plotted energy versus 
translation eigenvalue. One may clearly see from their figure that Eo{Ciq, 8) < 
Eq{Ciq, 7), giving one more example of violation of FOEL. We did not attempt 
to reproduce their picture here. Firstly, we were not able to numerically imple- 
ment the ring for 16 sites. But more importantly, the interested reader will do 
better to read their very beautiful paper, which is available online. See reference 

la- 

4 Relation to the "spectral curve" and violation 
of the string hypothesis in Bethe's ansatz 

These small calculations beg the question whether these violations of FOEL 
persist, so that there are always violations of FOEL between levels n— 1 and n 
a.t N = 2n for n > 2. We believe that is true. But they also beg the question 
whether this is a consequence of finite size corrections, or whether this has to do 
with the limiting "spectral" curves. This is supposedly answered by Sutherland 
[18], Dhar and Shastry T, and Barghecr, Beisert and Gromov [3^. 

Sutherland deduced that for sufficiently large spin deviations, asymptoti- 
cally, the density of the string state would have to become greater than 1 in 
order to satisfy the string hypothesis. But this is impossible because the roots 
must satisfy a trigonometric identity that has a minimum spacing of 1. There- 
fore, Sutherland claims that there is a value of p* G [0, 1/2], such that if one 
considers 

EoiCN,iN/2)^k), 
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for k > p^,N, and sufficiently large N, then one needs to calculate a different 
"spectral curve," than the vertical strings typically assumed to describe the 
Bethe roots in the so-called "string hypothesis." Sutherland claims that the 
determination of the spectral curve is an example of a classical Cauchy problem, 
and finds elliptic curves which give the correction to the string hypothesis. 

In the last displayed equation of his paper he gives the following parametric 
formula for the density of spin deviates and the energy density: 

d = (1/2) + a[{E{l/a)/K{l/a)) - l]/2 , 
e = 4,K{l/a)[2E{l/a) - (1 - {l/a^))K{l/a)] . 

Here the number of spin deviates is n = dN and the energy is £■ = e/N . The 
energy scaling is reasonable because the spectral gap is 7 = 0{1/N'^), and one 
expects the minimum energy of singlets to be on the order of 0{'jN) because 
one adds an amount of energy on the order of 0(7) for each "spin deviate". 
(See, for example, our forthcoming paper 117] for a rigorous derivation of some 
parts of the linear spin wave hypothesis which suggests this energy scaling.) 

Sutherland's formula involves K{k) and E{k), the complete elliptic integrals 
of first and second kind, respectively. Sutherland notes that to obtain d = 1/2 
means one takes the parameter a — >■ 00. Physically, the parameter a is related 
to the length of the spectral curve along which the Bethe roots accumulate. At 
half-filling of spin deviates this is supposed to diverge. But for small values of 
k, one has the asymptotic expressions: 

(2n)!! ) ' +■■■)' 
(2n-l)!!y fc2" \ 




^2 •4/ 3 V (2")!! / 2n-l ■"/ 

(4.8) 
Using this, one does recover d — 1/2 as a ^ 00, which means fc = 1/a — >■ 0. 
Also, one recovers Sutherland's formula that 

e -> TT^ as d -J- 1/2 , 

in the large N limit. 

The curve above refers only to the minimum energy among eigenvectors 
with total angular momentum n. According to Sutherland's surmise, this gives 
exactly the correct formula at half-filling of spin deviates, i.e., in the subspace 
with total spin equal to zero. In particular, this means that e —> tt'^ as d —> 1/2 
is correct within the context of the Bethe ansatz, using Sutherland's correction 
of the string hypothesis. Moreover, considering this curve for d < 1/2 but very 
close to 1/2, one sees that it is monotone decreasing. On the other hand, the 
first derivative is zero. Therefore, one needs to use a fine asymptotic analysis 
to see the persistence of the violation of the FOEL. 

Dhar and Shastry continued Sutherland's work, starting from the "low den- 
sity" regime. They found numerically, that the value of p^ where the string 
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solution breaks down is around 1/4. At low densities, Dhar and Shastry discov- 
ered that the Bethe ansatz equations are asymptotically satisfied by the Hermite 
points: roots of the Hermite functions. They also find that there is an exactly 
quadratic formula, as a function of the density of spin deviates, at low density. 
In fact, in Sutherland's notation, their formula is 

e = 4:n^d{l - d) . 

Interestingly, their arguments indicate that this formula is not merely asymp- 
totic in the limit rf — ;> 0. It should match with Sutherland's formula, if one 



includes the correct angular momentum curves generalizing (4.8) for d < 1/2. 
(We are very grateful to an anonymous referee for clarifying this point to us, 
which we had mis-reported on in an earlier version of this draft.) To us, Dhar 
and Shastry's results seem to be wholly reliable. 

Bargheer, Beisert and Gromov continued this line of investigation in [3]. 
They were motivated by new connections between models such as the spin-1/2 
Heisenberg model and string theory apparently using the AdS/CFT correspon- 
dence. But their analysis seems to be a very careful extension of Sutherland's 
orginal method. In particular they give many details of the so-called "spectral 
curves," that replace the vertical lines corresponding to the string solution, at 
high density. They also gave a concrete and simple numerical algorithm for 
searching for the Bethe ansatz equations. It amounts to fixing the spin deviate 
parameter fc, and varying the length of the ring, N, treating it as a parameter. 

For N ^ k, one is at "low density." Therefore, Dhar and Shastry's Hermite 
solution is approximately valid. This is the solution that Bargheer, Beisert and 
Gromov start with for a large value of N. Then they run a Newton method 
solver to recover the exact formula. At each step they decrease N, using the 
solutions of the Bethe ansatz equation at the last step as the starting point 
for a Newton scheme. As long as N is sufficiently large, one may decrease A'' 
in steps of size 1, proceeding along integers. But as N becomes smaller, this 
step size is too large to hope that the solution for A^ -|- 1 is a good starting 
point for Newton's method for A^. Then they recommend to slowly decrease 
N, treating it as a real, continuous parameter. Of course for non-integer values 
of N, the Bethe ansatz equations are not algebraic. But the Newton's method 
approach still seems to converge. In our implementations of their method (both 
in Mathematica where they originally ran it, and in Matlab) we found that the 
approach becomes chaotic as A^ approaches 2k from above. This seems to be 
part of the overall picture of a change of nature of this band of Bethe solutions 
as d approaches 1/2. 

The violations of FOEL for spin rings are tantalizingly close to these inter- 
esting results for the Bethe ansatz solution of the spin-1/2 Heisenberg model 
on rings. We hope to return to this question later, or to discover a definitive 
answer within the physics literature. 
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